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ABSTRACT 

The CODE microwave background temperature fluctuations and the abundance of local rich clusters 
of galaxies provide the two most powerful constraints on cosmological models. When all variants of the 
standard cold dark matter (CDM) model are subject to the combined constraint, the power spectrum 
of any model is fixed to ^ 10% accuracy in both the shape and overall amplitude. These constrained 
models are not expected to differ dramatically in their local large-scale structure properties. However, 
their evolutionary histories differ, resulting in dramatic differences towards high redshift. In particular, 
it should be true that any statistical measure that probes a rapidly diminishing tail of some distribution 
should provide a sensitive test at some sufficiently high redshift, when the objects in question are rare 
and hence in the tail. 

We examine in detail six standardized, CODE and cluster normalized GDM models with respect to a 
large set of independent observations. The observations include correlation function of rich clusters of 
galaxies, galaxy power spectrum, evolution of rich cluster abundance, gravitational leasing by moderate- 
to-high redshift clusters, Lya forest, damped Lya systems, high redshift galaxies, reionization of the 
universe and future GMB experiments. It seems that each of the independent observations examined is 
or potentially is capable of distinguishing between at least some of the models. The combined power of 
several or all of these observations is tremendous. Thus, we appear to be on the verge of being able to 
make dramatic tests of all models in the near future using a rapidly growing set of observations, mostly 
at moderate to high redshift. Gonsistency or inconsistency between different observed phenomena on 
different scales and/or at different epochs with respect to the models will have profound implications for 
theory of growth of cosmic structure. 

Subject headings: Gosmology: large-scale structure of Universe - cosmology: theory - galaxies: 
clustering - galaxies: formation - numerical method 


1. INTRODUCTION 

On one hand, the success of the COsmic microwave 
Background Explorer (GOBE; Smoot et al. 1992) can 
hardly be overstated; it revolutionized our way of quan¬ 
tifying cosmology. On the other hand, it did little justice 
to enthusiastic cosmologists: rather than helping us nail 
down to an ultimate model, it proded us to introduce more 
models (or variants) to meet with the new observational 
constraint. Post-GOBE era is filled with new ideas about 
models at the oxponoo of introducing more paramotoro to 


spectrum: the shape of the power spectrum (character¬ 
ized by the primordial power index on the large scales, n) 
and/or the composition of the universal matter including 
vacuum energy density (cosmological constant). It should 
be noted that variations of the latter also change the over¬ 
all shape of the power spectrum, quite separate from the 
primordial spectrum. The consequence of allowing such 
added freedom is the current state of a degenerate set 
of “viable” GOBE and cluster-normalized GDM models 
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for an excellent review). The flourish of new models is 
largely thanks to the general consensual requirement that 
any cosmological model match both the GOBE observation 
on large scales at an early epoch and local observations, 
primarily of distributions of galaxies. 

The newly introduced models that are deemed to be in¬ 
teresting are mostly based on the general success of the 
standard cold dark matter model. The variants arise by 
introducing one or two more parameters to the standard 
critical density GDM model with Zel’dovich primordial 


|Ostriker fc Steinhardt 1995|) and the open cold dark mat¬ 

ter model (Gott 1982; BucLct, Goldhaber, fc Turok 1995). 
Although noticeably being paid significantly less attention 
to, the primeval baryon isocurvature model (PBI, Peebles 
1987a,b; Gen, Ostriker, fc Peebles 1993) still appears con¬ 
tentious, especially with its virtue of early structure for¬ 
mation. 

After a general world model is chosen (i.e., flo, Aq, Hq 
and dark matter and baryonic compositions are chosen), 
the shape and the overall amplitude of its power spectrum 
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are yet to be fixed before the model becomes completely 
deterministic. COBE provides one fixing point at large 
scales (A ~ 1000/i“^Mpc). The abundance of local rich 
clusters of galaxies provides another fixing point at inter¬ 
mediate scales (A ~ 10/i“^Mpc). Both points are cur¬ 
rently determined to an accuracy of ~ 10%. Once the 
power spectrum is also fixed, the entire history of a model 
universe can, in principle, be computed. Sufficiently accu¬ 
rate numerical techniques with sufficiently realistic phys¬ 
ical modeling have become more available due to parallel 
advances in both computer hardware and software capa¬ 
bilities to allow for accurate predictions for at least some 
aspects or phenomena of the cosmic structure formation 
to be made, which can be compared to observations on 
various scales and at various epochs. Examples are Lya 
clouds (Cen et al. 1994; Zhang et al. 1995; Hernquist et 
al. 1996; MCOR), X-ray clusters (Kang et al. 1994; Bryan 
et al. 1994; Cen & Ostriker 1994a; Frenk et al. 1998) and 
gravitational leasing by clusters of galaxies and large-scale 
structure (Cen et al. 1994; Wambsganss et al. 1995). 

Parallel in time, rapid advances in observational capa¬ 
bilities in recent years, thanks in large part to the ground- 
based 10-meter class Keck telescopes as well as the Hub¬ 
ble Space Telescope (HST), make possible explore regimes 
which were previously unaccessible. For example, HST 
opens a brand new window to studying the formation and 
evolution of galaxies to a much higher redshift than before 
(e.g., Lowenthal et al. 1997) and quasar absorption sys¬ 
tems at some entirely new wavelength bands (e.g., Bahcall 
et al. 1996). Observations by Keck telescopes have made 
possible routinely record high quality quasar absorption 
spectra (e.g., Hu et al. 1995), providing new opportunities 
to conduct detailed studies of Lyman alpha forest lines 
and related metal line systems (e.g., Lu et al. 1997). Ob¬ 
servations of Lyman break galaxies using Keck telescopes 
have opened yet another new window to the universe at 
high redshift (Steidel et al. 1998, S98 henceforth). Deeper 
and larger surveys of X-ray clusters of galaxies as well as 
optically identified clusters of galaxies are pushing the red¬ 
shift of observed rich, luminous and hot clusters to beyond 
unity (e.g., Rosati et al. 1998). More excitingly, more and 
rapid observational advances are expected to take place 
in the near future in anticipation of many large projects 
and next generation telescopes convering a variety of wave¬ 
lengths from meterwave to y-ray. 

This paper attempts to interface between these two de¬ 
velopments. To set a standard for comparing different 
models and for comparing results between different work¬ 
ers in the field, we adopt the set of six CDM models pre¬ 
sented in Cen (1998; C98 henceforth) for detailed analyses. 
All the models are normalized to both COBE on large 
scales and the observed abundance of local rich clusters 
of galaxies, thus titled “viable”. The goal of this paper 
is to address the opportunities that new observations are 
and will soon be providing us with new ways to test these 
remaining, “viable” models. It is hoped that some theo¬ 
retical insight may be extracted to help highlight some po¬ 
tentially promising observations for (at least) the purpose 
of testing cosmological models. It is exciting that each 
of the independent observations examined here, mostly 
at moderate-to-high redshift, seems to be or potentially 
be capable of distinguishing between at least some of the 
models. The combined power of several of these obser¬ 


vations is tremendous, enabling the possibility of nailing 
down the ultimate cosmological model with high precision. 
It is perhaps more important that different observations at 
completely different scales and/or epochs should provide 
an overall consistency check of the general gravitational 
instability picture and/or Gaussianity of the primordial 
density fluctuations. 

The main results of the paper are organized as follows. 
The parameters of a suite of six CDM models that are 
normalized to COBE and local cluster abundance are dis¬ 
cussed in §2. §3 presents tests on the models with nine 
independent observations in order of ascending redshift 
in the subsections: the cluster correlation function (§3.1), 
the galaxy power spectrum from the Sloan Digital Sky 
Survey (SDSS) (§3.2), the evolution of rich cluster abun¬ 
dance (§3.3), the gravitational leasing by moderate-to-high 
redshift clusters of galaxies (§3.4), the Lya forest (§3.5), 
the damped Lya systems (§3.6), the high redshift galax¬ 
ies (§3.7), the reionization of the universe (§3.8) and the 
future CMB experiments (§3.9). Discussion of the con¬ 
straining/differentiating power and observational require¬ 
ment of each of the suite of observations is presented in 
§4. Conclusions are given in §5. 

2. SIX COBE AND CLUSTER NORMALIZED CDM MODELS 

The models are taken from C98 (Table 4 of that paper) 
and listed in Table 1. When normalized to both COBE on 
very large scales and the abundance of local rich clusters 
of galaxies, both the shape (n) and amplitude of the power 
spectrum (erg) of any model are fixed to about 10% accu¬ 
racy (C98). These six models are determined to this accu¬ 
racy. The tCDM model has Hq = 1 but with a tilt on the 
primordial power spectrum. The HCDM model has Hq = 1 
with 20% mass in massive neutrinos. Two neutrino species 
of equal mass are assumed. OCDM25 and OCDM40 are 
two open models and ACDM25 and ACDM40 are two low 
density but flat models with a cosmological constant. Ex¬ 
cept the open models, a tensor component is added when 
normalizing each model to COBE with tensor to scalar 
ratio T/S = 7(1 — n) (Liddle & Lyth 1992; Davis et al. 
1992; Crittenden et al. 1993; Stewart & Lyth 1993). Note 
that some tilt is present in all models listed, required to 
match both COBE and local cluster abundance. This set 
of six CDM models likely brackets all potentially interest¬ 
ing models. 

The baryon densities for OCDM25 and ACDM25 are 
computed using Vl^h? = 0.0125 (Walker et al. 1991), and 
for the remaining four models using = 0.0193 (Buries 
& Tytler 1997). These choices of Of, for the models are 
made to maximize the viability of each model with respect 
to the observed gas fraction in X-ray clusters of galax¬ 
ies (White et al. 1993b; Lubin et al. 1996; Danos & Pen 
1998, Pgas/Plot = (0.053 ± 0.004)/i“^/^). The power spec¬ 
trum transfer functions for all the models are computed 
using CMBFAST code developed by Seljak and Zaldar- 
riaga. The choice of the Hubble constant is made for 
each model such that each model is consistent with cur¬ 
rent measurements of the Hubble constant. It appears 
that H[j{obs) = 65± lOkm/s/Mpc can account for the dis¬ 
tribution of the current data from various measurements 
(see, e.g., Trimble 1997), except for those from Sunyaev- 
Zel’dovich observations (for a discussion of a reconcilia- 



3 


Table 1 

Six COBE and Cluster-Normalized CDM Models 


Model 

Ho 

n 



Ao 

fib 

crs 

tCDM 

55 

0.77 

0.936 

0.00 

0.0 

0.064 

0.55 

HCDM 

55 

0.88 

0.736 

0.20 

0.0 

0.064 

0.52 

OCDM25 

65 

1.47 

0.220 

0.00 

0.0 

0.030 

1.00 

OCDM40 

60 

1.15 

0.346 

0.00 

0.0 

0.054 

0.80 

ACDM25 

65 

1.10 

0.220 

0.00 

0.75 

0.030 

0.95 

ACDM40 

60 

0.96 

0.346 

0.00 

0.60 

0.054 

0.80 


Figure 1 



Fig. 1.— Cluster-cluster two-point correlation length computed using GPM as a function of Qq for clusters with mean separation of 
55/i~^Mpc (light dashed curves) and 31/i~^Mpc (light solid curves), respectively. For all models we use the CDM transfer function from 
BBKS with F = Qoh = 0.23 and n = 1.0. All models are normalized to clusters using equation 36 of C98. Each set of four curves consists of 
two upper curves using 2cr upper limit of as by equation 36 of C98 with Aq = 0 and with Qq + Aq = 1, and two lower curves using 2a upper 
limit of as with Aq = 0 and with Hq + Aq = 1. Also shown as two heavy dashed horizontal line indicate the 2a upper and lower bounds from 
observations by Peacock & West (1992) of the correlation length of clusters of mean separation of 55/i~^Mpc. The two heavy solid horizontal 
lines are the 2a upper and lower bounds from observations by Dalton et al. (1994) of the correlation length of clusters of mean separation of 
31/i“^Mpc. 
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tion of this difference, see C98). Another consideration 
is that the age constraint from latest globular cluster ob¬ 
servations/interpretations (c.f., Salaris, Degl’Innocenti, & 
Weiss 1997) is not violated. 

3. TESTING COM MODELS AT MODERATE TO HIGH 
REDSHIET 

While COBE and local cluster abundance combine to 
provide an unprecedented accuracy of fixing models in 
terms of limiting the parameter space, some degeneracy 
remains. It is therefore essential to examine other, inde¬ 
pendent observations to help further tighten the parameter 
space. We will examine various observations in this sec¬ 
tion using the six models listed in Table 1 as a testbed. We 
consider nine independent observations in order of ascend¬ 
ing redshift: the cluster correlation function at low red- 
shift (§3.1), the galaxy power spectrum from SDSS (§3.2), 
and evolution of rich cluster abundance (§3.3), the grav¬ 
itational leasing by moderate-to-high redshift clusters of 
galaxies (§3.4), the Lya forest (§3.5), the damped Lya 
systems (§3.6), the high redshift galaxies (§3.7), the reion¬ 
ization of the universe (§3.8) and the future CMB experi¬ 
ments (§3.9). 

3.1. Constraints by the Cluster Correlation Function at 

z = 0 

Before examining the correlation functions of the six 
models, it is worthwhile to understand factors that deter¬ 
mine the cluster correlation function. Figure 1 shows the 
cluster-cluster two-point correlation length (i.e., the scale 
at which the correlation function is unity), computed using 
the Gaussian Peak Method (GPM, G98), as a function of 
flo for clusters with mean separation of 55/i“^Mpc (light 
dashed curves) and 31/i“^Mpc (light solid curves), respec¬ 
tively. For all models shown here, we use P = 0.23 with 
the GDM transfer function from BBKS and normalize to 
the local cluster abundance (equation 36 of G98) but not 
necessarily to GOBE. From Figure 1 it immediately be¬ 
comes clear that the cluster correlation function almost 
does not depend on Aq at all, only weakly depends on Oq 
and somewhat moderately depends on (jg: dro/dAo « 0, 
dro/dflo « 2 — 4, and dro/das « 6 — 10. Therefore, given 
the uncertainties in fixing as of a model, uncertainties in 
the shape of the power spectrum and uncertainties in the 
observed correlation length of clusters, it is not possible 
to constrain fig using correlation functions of clusters of 
galaxies, unless models with different fig also substantially 
differ in Pk on the relevant large scales (ii0 — 300h~^Mpc). 

Fortunately, the power spectra of various models with 
different flg/Ag are indeed quite different (Table 1), as 
required to fit both COBE and local rich cluster abun¬ 
dance. Figure 2 shows the power spectra (in real space) 
for the six models in Table 1. It is seen that normaliz¬ 
ing to both COBE and clusters results in a wide disparity 
in shape and amplitude of the power spectra of the mod¬ 
els in a wide range of scales. The correlation length of 
clusters of galaxies in the six models are then shown in 
Figure 3 in a somewhat different manner, where the cor¬ 
relation length in each model is plotted against the mean 
separation of the clusters in question. Note that the mean 
separation of a set of clusters is the mean separation of a 
set of the richest (most massive) clusters in a model; clus¬ 
ters are rank-ordered in mass and selected in descending 


order in mass. We see that there exists significant dif¬ 
ferences between different models shown here. Since we 
have explicitly shown earlier that the dependence of the 
correlation length on fig and Ag is rather weak, the differ¬ 
ences in the correlation lengths between the models shown 
here chiefly attribute to the differences in the power on 
large scales between the models. The trend is clear and 
physically understandable: the more power a model has 
on large scales A > 100h“^Mpc, the stronger dependence 
of the correlation length on the richness of the clusters 
(see Figure 2 to compare the power spectra of the models) 
is displayed, consistent with previous findings (Bahcall & 
Gen 1992; Croft & Efstathiou 1994). 

It is evident that the correlation length of clusters of 
galaxies sets a useful constraint on models. If the measure¬ 
ment accuracy can reach < lh“^Mpc, one can potentially 
strongly differentiate between models. It is also clear that 
it will be more useful to reduce the observational error- 
bars of the correlation length of very rich clusters (e.g., 
d > 90/i“^Mpc) for the purpose of testing models, since 
the models differ most there. Stated somewhat more quan¬ 
titatively in a hypothetic sense, were the 2a errorbar of the 
solid dot in Figure 3 smaller than 5/i“^Mpc, all the mod¬ 
els considered here are ruled out at a very high confidence 
level (~ Aa). 

Consistent with other observations (e.g., Maddox et al. 
1990) and previous results of clustering analyses of galaxy 
clusters (e.g., Bahcall & Gen 1992; Croft & Efstathiou 
1994), the standard Og = 1 GDM model even with a sig¬ 
nificant tilt of the power spectrum {n = 0.77, tCDM) is in 
disaccord with currently available observations due to the 
lack of large-scale power. HCDM model and all other, low 
density models are consistent with available observations 
of cluster correlations. 

3.2. Calaxy Power Spectrum 

The Sloan Digital Sky Survey now underway (Knapp, 
Lupton, & Strauss 1996) will provide a tremendously 
large database for a variety of objects. Realizing its full 
value will perhaps require the work by generations of as¬ 
tronomers. One obviously valuable piece of information 
that can be extracted immediately upon the availability 
of the distribution of galaxies from SDSS, up to very large 
scales previously unreached in full three dimension. Figure 
4 shows the power spectra of galaxies for the six models 
(Table 1) in redshift space at z = 0. The shown galaxy 
power spectrum in redshift space for each model is related 
to the linear matter power spectrum in real space by the 
following relation: 

Pgal,red{k) = (1 -f ^ ^)'^8 Pmass,real{k)- (1) 

In the above relation the first term on the right hand side 
is due to linear peculiar velocity distortion (Kaiser 1987) 
and the second term on the right hand side reflects the 
(linear) bias of galaxies over mass in order for each model 
to be consistent with the observed galaxy density fluctu¬ 
ations at z=0. For tGDM, we also show the errorbars 
{2a) that SDSS can achieve, adopted from Goldberg & 
Strauss (1998). Errorbars for other models are the same 
but not shown for readability of the figure. It is clear 
that the SDSS galaxy power spectrum in a large range 
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Fig. 2.— shows the power spectra (in real space) for the six models in Table 1. 


Figure 3 
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Fig. 3.— shows the correlation length against the mean separation of the clusters in the six models (Table 1). Shown as symbols are ob¬ 
servations from various sources, where the filled dot is for R> 2 Abell clusters from Bahcall Soneira (1983) with separation d = 94h~^Mpc 
and correlation length vq = 42 di 20h“^Mpc, the open circle is for R > 1 Abell clusters from Peacock Sz West (1992) with separation 
d = 51.8h“^Mpc and correlation length ro = 21.3 ib 2.6/i~^Mpc, the solid square is for Edinburgh-Durham Southern Galaxy Catalogue 
(EDSGC) from Nichol et al. (1992) with separation d = 46.4h~^Mpc and correlation length ro = 16.4di8.0h~^Mpc, and the two open squares 
are for APM Galaxy survey from Dalton et al. (1992) with separation d = 34.7h~^Mpc and correlation length ro = 12.9 di 2.8h“^Mpc and 
from Dalton et al. (1994) with separation d = 31.0/i“^Mpc and correlation length ro = 14.3 di 2.4/i~^Mpc. All the observational errorbars 
are 2a. 
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Figure 4 
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Fig. 4.— The power spectra of galaxies in redshift space for the six models (Table 1) at 2 = 0. See text for how they are obtained. For 
tCDM, we also show the errorbars (2cr) that SDSS are thought to achieve, adapted from Goldberg & Strauss (1998). 


in k from 0.01 ft./Mpc to 0.5 /i/Mpc will provide some 
discriminatory tests of the models. Taken at face value, 
it seems that the errorbars are smallest at the highest k, 
where the data points may provide the most discriminat¬ 
ing test of the models. But we caution that nonlinear 
effects may wash out some of the differences between the 
models at the high k end (e.g., Gramann, Cen & Bah- 
call 1993). We should only conservatively estimate the 
discriminating power of the galaxy power spectrum from 
SDSS at fc < 0.2 /i/Mpc, where the nonlinear effects are 
relatively mild. A close examination of Figure 4 reveals 
that there is not a pair of models that have the same 
power at fc < 0.2 /i/Mpc within the 2a errorbars provided 
by SDSS. Therefore, SDSS should be able to test all these 
models at a confidence level > 2a. The purpose of this 
investigation is to illustrate the potential of SDSS galaxy 
power spectrum in placing strong constraints on models. 
Quantitative comparisons will have to be made between 
large-scale, fully nonlinear cosmological galaxy formation 
simulations (Cen & Ostriker 1992b,1993b,1998) and obser¬ 
vations. 

3.3. Evolution of Cluster Mass Function 

Let us now study the evolution of the cluster mass func¬ 
tions in the six models (Table 1), computed using GPM 
(C98). Figure 5 displays the mass functions for the models 
at four redshifts, 2 ; = 0.0, 0.3, 0.5,1.0. The cluster mass is 
defined as the mass in a sphere of comoving Abell radius 
IjA = 1.5 ft.“^Mpc) at all redshifts. Again, all models are 
normalized to the observed mass function at redshift zero 
as well as COBE. As one would have naively expected, 
the evolution is strongest in the Dq = 1 models (5a,b) and 


weakest in the open, lowest density model (5c). It must be 
stressed, however, as we will show shortly, that the phys¬ 
ical reason for this trend may come as a surprise. The 
underlying cause is as, not flo- The HCDM (5b) shows 
a comparable rate of evolution to tCDM (5a) because the 
reduction of power on smaller scales k > 0.12h/Mpc is 
roughly compensated by a comparable increase of power 
on larger scales k < 0.12h/Mpc in HCDM compared to 
tCDM. Other models show intermediate rates of evolution, 
in a comprehensible order. Figure 6 is similar to Figure 
5 with the only difference being that the cluster mass is 
now defined as the mass in a sphere of proper Abell radius 
ijA = 1.5 /i“^Mpc) at all redshifts. The trend between 
models is similar in Figure 6 to that in Figure 5. But the 
evolutionary rates in all models are dramatically smaller 
in Figure 6 than in Figure 5. This says that the mass of 
clusters within a physical radius evolves much more slowly 
than the mass of clusters within a comoving radius. This 
is trivial, since a fixed physical radius corresponds to a 
larger comoving radius at a higher redshift thus its sphere 
contains more mass than its counterpart with the comov¬ 
ing radius. A more subtle difference between the models 
is that the evolution of the mass function for the low den¬ 
sity open model (OCDM25) to redshift unity (the bottom 
curve in Figure 6c) is much smaller than that of any other 
model. In other words, the mass of clusters within a physi¬ 
cal radius in OCDM25 evolves only very slowly. An equiv¬ 
alent, conventional but perhaps conceptually simpler way 
to put this is that the clusters in OCDM25 virialize early 
at 2 : ~ 1 and later mass accretion onto the formed clusters 
through the surface of the indicated radius is modest. 

We would summarize the trends for clarity. Clus- 
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Figure 5 
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Fig. 5.— It displays the mass functions for the models at four redshifts, z = 0.0, 0.3, 0.5,1.0. The cluster mass is defined as the mass in 
a sphere of comoving Abell radius = 1.5 h“^Mpc) at all redshifts. All models are normalized to the observed mass function at redshift 
zero. 
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ter mass function is a monotonic function in time: the 
abundance of clusters at a fixed mass (at least in the 
range of mass considered here from 4.8 x M© to 

1.2 X decreases with increasing redshift. For 

all COBE and cluster (at z = 0) normalized models, the 
lower is Hq, the weaker is the evolution; for two models 
with a same Hq, an addition of a cosmological constant 
enhances the evolution, especially beyond a certain high 
redshift when the values of in the two models (with 
and without a cosmological constant) start to deviate sig¬ 
nificantly from one another. However, it is worth noting 
that the faster evolution in the low density models than in 
the high density models is somewhat coincidental (due to 
(Tg —Hq relation required to fit local rich cluster abundance, 
equation 36 of C98) and entirely due to the normalization 
on (Tg, not due to Hq, as will be shown below. 

Presenting the above results shown in Figures (5,6) in 
a slightly different form, panel (a) of Figure 7 shows the 
evolution of the number of clusters with mass of > 5.6 x 
M© within a fixed proper radius of 1.5/i“^Mpc, 
and panel (b) of clusters with mass of > 7.9 x M©. 

The quantity in the ordinate is the number of clusters per 
unit redshift per unit steradian. Since the number of clus¬ 
ters per unit redshift per unit sold angle is plotted, one has 
to take account the different geometries in different model 
universes (the comoving volume element per steradian per 
unit redshift about 0 is DA{z)'^drcm{,z)/dz, where Da is 
the comoving angular diameter distance to redshift z and 
fcmiz) is the comoving distance to redshift z; see equations 
6, 8 below for definitions of Da and rcm)- It is seen that 
the less massive, low redshift clusters shown in panel (a) 
has much less discriminating power than the more mas¬ 
sive, higher redshift clusters shown in panel (b). While 
the observed clusters at z = 0.18 — 0.35 (panel a) are not 
able to differentiate between any models at a significant 
confidence level, the observed clusters at z = 0.35 — 0.55 
and z = 0.55 — 0.85 (panel b) offer a clean way to separate 
the Hq = 1 models from the rest. It seems quite secure 
to conclude that both Hq = 1 tCDM and HCDM models 
are ruled out at a high confidence level (> 3a), whereas 
the current observations are not able to differentiate be¬ 
tween the four low density models. This conclusion is in 
broad agreement with the conclusions drawn by Henry & 
Arnaud (1991), Bahcall, Fan & Cen (1997) and Carlberg 
et al. (1997). Reichart et al. (1998) claim that the X-ray 
cluster luminosity function from Einstein Medium Sensi¬ 
tivity Survey are consistent with Hq = 1. It is unclear at 
this time what causes the difference. However, Reichart 
et al. (1998) claim that their conclusion is independent of 
(Tg (i.e., the amplitude of the power spectrum), which is 
in direct conflict with the demonstration that CTg is the 
primary factor that determines the redshift evolution of 
massive clusters, as shown below. 

To quantify the power of high redshift X-ray clusters 
for differentiating between models, we tabulate in Table 2 
the expected number of clusters in each of the six mod¬ 
els at redshifts greater than 0.5, 1.0 and 2.0, respectively, 
for clusters with temperatures greater than 0.5keV and 
7.0keV, respectively. The numbers are obtained assuming 
that a survey covering tt steradian (i.e., one quarter of the 
sky) is available. Next generation X-ray telescopes should 
be able to detect 7keV clusters at z = 1 without much 
difficulty, if clusters at high redshift are not significantly 


less luminous than their local counterparts at a fixed tem¬ 
perature (e.g., Henry 1997). It is clear from Table 2 that 
Ho = 1 models are easily differentiated from other models. 
The subtlety lies in low density models and some degener¬ 
acy exists between open models and flat models with cos¬ 
mological constants. For example, OCDM40 (Hg = 4.0, 
Ao = 0.0) and ACDM25 (Hg = 0.25, Ag = 0.75) give 
comparable number of bright X-ray clusters at z < 1.0; 
but the difference is gradually larger at higher redshift (a 
factor about four at z = 2). Aside from this degeneracy, 
the database from a survey of the indicated sky coverage 
should be able to determine Hg with an uncertainty of 
AHg = 0.02 (Ict), if Hg -I- Ag = 1 [obtained by assuming 
statistical (Poisson) errors for the clusters with kT = 7keV 
at z > 1 of the last two models in Table 2, which exhibit 
the smallest difference and thus are most demanding]. If 
the universe is open, one will be able to determine Hg with 
an accuracy of AHg = 0.008 (by interpolating between 
OCDM25 and OCDM40). These accuracies are compara¬ 
ble to those which will be determined from the next gener¬ 
ation microwave background fluctuation experiments (e.g., 
Zaldarriaga, Spergel, & Seljak 1997), thus provide a very 
powerful test and also a very important consistency cross¬ 
check. 

So far we have shown that there are dramatic differ¬ 
ences in the number density of massive clusters at high 
redshift between the considered models, when all are nor¬ 
malized to give the same cluster density at z = 0. In 
other words, the rate of evolution for cluster abundance is 
strongly model-dependent. It is useful to understand what 
cause such dramatic differences in the evolution rates be¬ 
tween the different models. For this let us first define an 
evolution rate of cluster (with mass > M) number density 
as a function of as as 

( 2 ) 

dz 

This is the change of cluster abundance in logarithmic 
scale per unit redshift. This definition is similar to that 
defined in Fan, Bahcall, & Cen (1998; equation 8) but with 
the difference that the current £ is a function of as and 
M, not a function of z and M. In this sense the current 
definition of rate of evolution is local in time; redshift z in 
equation 2 serves as a ruler or a unit. Figure 8a shows the 
evolution rate of the number density of clusters with mass 
> 5.6 X M©, defined within a comoving Abell ra¬ 

dius (panel a) and within a physical Abell radius (panel b), 
respectively. The first thing to notice is that, not surpris¬ 
ingly, the evolution rate is lower for clusters defined within 
a physical Abell radius than within a comoving Abell ra¬ 
dius. Secondly, the differences between models at a fixed 
as are smaller than the differences of any single model at 
different as’s. For example, for OCDM25 S changes from 
1.4 to 6.8 for as from 1.0 to 0.5, while the differences in S 
between tCDM and OCDM25 is 0.6 at as = 1-0 and 2.0 
at (Tg = 0.5. However, it is interesting to note that, at 
a fixed as, the evolution rates in the low Hg models are 
higher, not lower, than the rates in the Hg = 1 models for 
clusters defined within the comoving Abell radius (pan¬ 
els a). The situation is altered for clusters defined within 
the physical Abell radius: at a fixed as, the evolution rate 
is lowest in OCDM25. In any case these differences are 
secondary compared to the difference of £ between differ- 
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Figure 7 


1000 


1000 



100 


10 


1 


0.1 



Fig. 7. — Panel (a) shows the number of clusters with mass > 5.6 X Mq within a fixed proper radius of 1.5/i“^Mpc, and panel (b) 

shows the evolution of the number of clusters with mass > 7.9 X Mq within a fixed proper radius of 1.5/i“^Mpc. The quantity in 

the vertical axis is the number of clusters per unit redshift per steradian. Shown as symbols are the observations with 2cr errorbars (Carlberg 
et ah 1997, C97; Luppino & Gioia 1995, LG95). 


Table 2 

Number of High Redshift X-ray Clusters®- in Six Models 


Model 

(5.0,7.0)keV 

0 > 0.5 

(5.0,7.0)keV 
z > 1.0 

(5.0,7.0)keV 
z > 2.0 

tCDM 

(71,2.2) 

(1.3,0.01) 

(0.0,0.0) 

HCDM 

(24,1.3) 

(0.19,0.0028) 

(0.0,0.0) 

OCDM25 

(14902,1374) 

(10591,798) 

(3232,139) 

OCDM40 

(3323,273) 

(1314,70) 

(75,1.4) 

ACDM25 

(3217,369) 

(1005,72) 

(16,0.30) 

ACDM40 

(1270,112) 

(203,9.1) 

(0.56,0.0040) 


-for a sky coverage of tt steradian 
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Fig. 8 .— ( 8 a) shows the evolution rate of cluster number density as a function of the amplitude of the power spectrum (crsl for the six 
models (Table 1). Panels (a,b) of Figure ( 8 a) show the evolution rate of number density for clusters whose mass (=5.6 x lO^‘^/i~^M 0 ) is 
defined within a comoving Abell radius and within a physical Abell radius, respectively. ( 8 b) shows the dependence of the evolution rate on 
power index n by varying n for tCDM and ACDM25. 
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ent (7s within a same model. Next, we select two models 
(tCDM and ACDM25) and let the power index n vary 
within the allowable ranges to test the dependence of £ 
on n. The results are shown in Figure 8b, from which we 
can easily draw the conclusion that the dependence of the 
evolution rate on n is very weak. 

We conclude that the dramatic differences in the evolu¬ 
tion rates between the different models are almost entirely 
caused by the differences in erg between the models, which 
results from the requirement that each model has to fit 
the cluster abundance at z = 0. The apparent dependence 
of the evolution rates on flo/Ao is dictated by the strong 
dependence of cluster abundance on erg through the erg — flo 
relation (equation 36 of C98). This conclusion is in good 
agreement with that of Fan et al. (1998), who use the 
Press-Schechter method. 


3.4. Gravitational Tensing by Clusters of Galaxies 

Clusters of galaxies are so massive that they ubiqui¬ 
tously cause noticeable distortions of images of back ground 
sources even in their outskirts (~ 1 /t~^Mpc) (iTyson 


Valdes, fc Wenk 1990; Miralda-Escude 1991; Kaiser & 


Squires 1992). The central regions of rich clusters are often 


capable of producing more spectacular events such as large 
splitting, multiple images (of background quasars) or giant 
arcs (of background galaxies) (e.g., Le Fevre et al. 1994; 
Edge et al. 1994; Small et al. 1996). These extreme grav¬ 
itational leasing events are rare because the cross section 
for such events is small; usually only the central region of 
size of order several kiloparsecs in a very rich cluster is ca¬ 
pable of producing giant arcs or multiple splittings of size 
of arcseconds and larger, compared to the usual notion of 
the size of a cluster (~ IMpc) or even the size of the X-ray 
core (~ lOOkpc) of a cluster. Consequently, the frequency 
of such events is extremely sensitive to abundance of a 
very small population of very rich clusters at intermediate 
redshifts (^ 0.3 — 2.0), where leasing is optimal. 

Narayan & White (1998) compare strong leasing events 
in various models (pre-COBE) with observations combin¬ 
ing Press-Schechter method with the singular isothermal 
sphere assumption. Cen et al. (1994) and Wambsganss 
et al. (1995) use detailed N-body simulations to compute 
the leasing statistics, focusing on intermediate splittings (a 
few arcseconds). Here we compute the frequency of very 
large splitting or giant arc events (~ areminute). We 
approximate each cluster as a singular isothermal sphere 
(SIS; Turner, Ostriker & Gott 1984, TOG henceforth) with 
the one-dimensional velocity dispersion determined by the 
mass within the comoving Abell radius: 


2 GAf^(l -|- zl) 

tTn — ■ 




(3) 


where G is the gravitational constant, and Zl is the red- 
shift of the (lens) cluster in question. For a SIS the bending 
angle, 9q, conveniently independent of impact parameter, 
is (TOG) 

eo = M—?, (4) 

c 

where c is the speed of light. The size of the arc due to 
such a lens will be approximated as the diameter of the 
Einstein ring as: 


29oDls 
Ds ’ 


where D^s and Ds are the (comoving) angular diameter 
distance between the lens (cluster) and the source, and 
between observer and the source, respectively. The (co¬ 
moving) angular diameter distances are dependent upon 
the geometry of a world model [D^s = D{xs — Xl) and 
Dsixs)] as: 

{ hoX for flat universe (Hq + tVo = I) 

qq sin(x) for close universe (Hq -I- Aq > I) 
qo sinh(x) for open universe (Hq -I- Aq < l) 

( 6 ) 

Here, uq is a dimensional scale factor defined below; x is 
the dimensionless comoving distance to the object in ques¬ 
tion at z, which can be expressed for an arbitrary universe 
(Refsdal, Stabell & de Lange 1967; Gott, Park & Lee 1989) 
as 

X{z) = \/|3(To - qo - 1| [ Z{u)du, (7) 

do 

where Z(u) = [l+2(qo + l)u+(qo + l + 3ao)u^+2aou^]~^^^. 
Here cto = ^o/2 and qg is the deceleration parameter, 
which is qo = Ho/2 for Aq = 0 and qo = 3Ho/2 — 1 
for Hq -I- Aq = 1. The dimensional scale factor oq is 
oo = - , , where Hq is the Hubble constant. 

IS'^o—90 —1| 

The dimensional comoving distance of an object at z is 

rcm(z) = aox(z)- (8) 

Combining equations 3,4,5, we obtain an expression for 
the required cluster mass Ma in order to produce an arc 
of size [3 as 


MA{f3) 


(dc^VA Ds 
47rG(l -I- zl) Dls 


(9) 


We can then compute the optical depth per unit redshift 
for forming arcs of size > /?, as a function of redshift: 


dr 2 
dz dz ^ 


( 10 ) 


where is the number of clusters with mass > Ma{I3) 
(equation 9) per unit steradian per unit redshift, computed 
using GPM (C98). 

Figures (9a,b) show the (differential) gravitational leas¬ 
ing probability (i.e, differential optical depth) per unit red¬ 
shift as a function of redshift for arcs or splittings greater 
than 30 and 60 for sources at Zg = 1 and Zg = 3, respec¬ 
tively. In accord with the evolution of cluster abundance 
as shown in Figures (5,6,7), the Hq = I models (tCDM, 
HCDM) have the lowest optical depths due to the rapid 
decrease in the number of rich cluster towards high red¬ 
shift. The lowest Hq models (OCDM25, ACDM25) have 
the highest optical depths. However, the difference be¬ 
tween OCDM25 (open Hq = 0.25) and ACDM25 (flat 
Ho = 0.25) (similarly between OCDM40 and ACDM40) 
is more subtle: for sources at Zg = I ACDM25 has a 
somewhat higher optical depth than OCDM25 because the 
prolonged path length due to the cosmological constant in 
ACDM25 plays a significant role while the evolution rate 
of cluster abundance in the two models are comparable 
out to redshift ^ 0.5; but for sources at Zg = 3 OCDM25 
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Fig. 9.— The differential gravitational lensing probability (i.e., differential optical depth) per unit redshift as a function of redshift for arcs 
or splitting angles of multiple image events with sizes greater than 30 and 60 for sources at z = 1 (9a) and 2 = 3 (9b), respectively. 
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has a significantly higher optical depth than ACDM25 be¬ 
cause the much lower rate of cluster abundance evolution 
in OCDM25 dominates over other factors. 

The assumption of SIS is admittedly approximate, but 
used because of its analytical simplicity and because of its 
astrophysical relevance as used in many statistical studies 
(e.g., Tyson 1983; TOG; Narayan & White 1988; Fukugita 
& Turner 1991; Mao 1991). A perhaps more realistic pro¬ 
file is NFW profile (Navarro et al. 1996), which goes as 
r~^ in the center and as r~^ in the outskirt. Which pro¬ 
file, NFW or SIS, is more efficient for lensing depends on 
the core radius of clusters. If the core radius is sufficiently 
small, the NFW profile could be more efficient than a SIS 
profile, for a given cluster mass within a larger, given ra¬ 
dius (Maoz et al. 1997). On the other hand, if the SIS 
assumption overestimates the strong lensing probability, 
this would diminish the already very small lensing opti¬ 
cal depth in the high Oq models further (the reductions 
in model with higher optical depths are smaller), widen¬ 
ing the difference between models. Given the simplicity of 
our treatment, we choose to focus on the gross frequency 
of large gravitational arcs or large splitting multiple im¬ 
ages, which are relatively robust, not on the details such 
as length-to-width ratio of arcs. The latter may sensitively 
depend on the small scale structures in the clusters such 
as core radius, inner density slope and substructure (see 
Bartelmann et al. 1997 for an excellent study on arcs of 
smaller sizes). The lack of or inadequate treatment of 
galaxy formation and inaccuracy in the small scale density 
distributions (even with very large, pure N-body calcula¬ 
tions currently available) do not yet permit a sufficiently 
accurate characterization of lens potentials on small scales. 

Clearly, the gravitational lensing optical depths of large 
arcs (or large splitting images) for different models deviate 
rapidly from one another towards high redshift. We argue 
that the errors introduced by our approximate method are 
likely to be substantially smaller than the difference be¬ 
tween the models at moderate-to-high redshift. Detection 
of even a small sample (or non-detection) of such large 
lenses at redshift >1.0 and/or a large sample of high red¬ 
shift sources with such large arcs will provide a useful con¬ 
straint on the models. Systematic surveys for large-image 
separation lensed quasars (Maoz et al. 1997) or arcs are 
beginning to provide useful statistical samples, and larger 
data sets such as from SDSS will be extremely powerful 
in constraining cosmological models. With the data set of 
10® quasars from SDSS (assuming to be at = 3 for this 
simple illustration), we expect to see (inferred from panel 
a of Figure 9b) (1900 ± 44, 710 ± 27, 570 ± 24, 330 ± 18, 
lOOilO, 39±6) multiple-image quasars with image separa¬ 
tion > 30” in (OCDM25, OCDM40, ACDM25, ACDM40, 
tCDM, HCDM), respectively (where the uncertainties are 
Icr Poissonian). This would allow an unambiguous test of 
these six models. More likely, it would put a constraint on 
Do with ADq ~ 0.01 if the universe is open or ADq ~ 0.03 
if the universe is flat. A more stringent test can be ob¬ 
tained if one also has the redshift distribution of the lenses. 
One 


3.5. Lya Forest 

The abundant low column density Lya forest lines in 
quasar absorption lines provide a precise and unbiased 
tracing of the distributions of neutral hydrogen, in a va¬ 
riety of environments, up to the highest redshift where 
quasars are observed to exist (e.g., Rauch et al. 1992; 
Schneider et al. 1993; Hu et al. 1995; Tytler et al. 1995; 
Lanzetta et al. 1995a; Bahcall et al. 1996). Until very 
recently, models of quasar absorption systems were post 
hoc. Successfully putting the formation and evolution of 
Lya clouds into a cosmological context, computed ab initio 
from standard cosmological initial conditions (Gen et al. 
1994; Zhang et al. 1995; Hernquist et al. 1996; Miralda- 
Escude et al. 1996, MCOR hereafter; Bond & Wadsley 
1997), has set forth a brand new way to study these ob¬ 
jects. The surprisingly good agreement between simula¬ 
tions and observations in terms of matching the distribu¬ 
tions of line width parameter b and column density, among 
other quantities, clearly demonstrates that the basic cos¬ 
mological origin for Lya clouds is no longer questionable. 

It is to distinguish between competing cosmological 
models that interests us most here. We attempt to check 
illustratively if Lya forest may be used to test cosmologi¬ 
cal models. We will only examine models at z = 2, where 
good observational data are available. Figure 10 shows the 
power spectra (panel a) and the rms density fluctuations 
(panel b), both in real space and computed using linear 
theory, at the relevant scales {k ~ 1.0 — 30.0h/Mpc and 
H{z)rTH ~ 10.0 — 300.0km/s) at z = 2. In order for the 
models to be shown and compared on a common, relevant 
physical scale, the Hubble velocity instead of proper (or co¬ 
moving) physical distance is shown in the x-axis in panel 
(b). Also shown as a cross-shaded region is the observa¬ 
tional constraint adapted from Rauch et al. (1997). The 
adopted observational constraint is based on the range of 
the fluctuation amplitudes on the relevant scales spanned 
by the two models (SCDM and ACDM; Figure 8 of Rauch 
et al. 1997) examined by Rauch et al. (1997), which fit the 
observations well and provide a useful indication in the 
needed power on the relevant scales. The plotted vertical 
range of the shaded box is three times the range between 
SGDM and ACDM in Rauch et al. (1997), centered on the 
mean amplitude (at each x-axis value) averaged over the 
SCDM and ACDM models. 

The disparity of powers between the models on these 
scales relevant to the formation of Lya clouds is rather 
striking. The very slow growth of perturbations at low 
density open models is clearly manifested: OCDM25 (open 
Do = 0.25) has distinctly more power on all the shown 
scales than other models, while HCDM has power at 
k = 30/i/Mpc lower than that of OCDM25 by almost 
two orders of magnitude. These differences should result 
in distinctly different populations of Lya forest lines and 
flux distributions. We expect that the overall flux distri¬ 
butions (Rauch et al. 1997) of the models will be different, 
even with the liberty to adjust the ionization rate to fit 


the observed mean flux dec rement (Zuo & Lu 1993; Press 


XTL nm/rTn” %bicki, fc Schneider 1993]) . In particular, the low flux end 


at z > 2 in (OCDM25, OCDM40, ACDM25, ACDM40, 

tCDM, HCDM), respectively, that produce arcs or split¬ 
tings of size > 30”. We see that the differences between 
models in this measure are even more dramatic. 


of the 1-point flux distribution will be dramatically differ¬ 
ent between the models (Cen 1997a). Since the amount of 
observational data available (and accumulating) on Lya 
forest lines are very large, measures based on flux distri- 
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Fig. 10.— The power spectra (panel a) and the rms density fluctuations (panel b), both in real space, at the relevant scales 
(fc ~ 1.0 — 30.0h/Mpc and H[z)rTH 10.0 — 300.0km/s) at 2 : = 2, computed using linear theory. Also shown as a cross-shaded region 
is the observational constraint adapted from Rauch et al. (1997). In panel (b) the x-axis is in km/s, obtained by multiplying the physical 
scale {rj'fi) by the Hubble constant at 2 = 2, H{z). 


buttons (MCOR; Cen 1997a; Miralda-Escude et al. 1997) 
may enable tests of models with high statistical signifi¬ 
cance. Preliminary comparisons between models and ob¬ 
servations using the one-point flux distribution probability 
function (Rauch et al. 1997) seem to most favor the two 
A based models (ACDM25 and ACDM40), but disfavor 
HCDM and OCDM25 in the opposite sense that HCDM 
has too little power and OCDM25 has too much power 
on the relevant scales. tCDM and OCDM40 appear to be 
marginally consistent. 

Definitive theoretical predictions of all the models and 
comparisons between models and observations await de¬ 
tailed modeling using cosmological hydrodynamic simula¬ 
tions suitable to study the formation of Lya clouds (Cen 
et al. 1994; Zhang et al. 1995). 

3.6. Damped Lya Systems 

The high end of the Lya clouds is the damped Lya sys¬ 
tems, which are thought to be progenitors of large disk 
galaxies (Wolfe 1995). Damped Lya systems are much 
rarer than the common Lya forest lines by approximately 
two orders of magnitude, analogous to relative abundances 
between galaxies and clusters of galaxies at present time. 
In this sense damped Lya systems probe the tail of the 
density distribution on ~ lh“^Mpc scale at high redshift, 
whereas present clusters of galaxies probe the tail of the 
density distribution on ~ 10h“^Mpc scale. 

Prochaska & Wolfe (1997) And that the only model 
among tested that matches the observed kinematics of 
damped Lya systems is a thick fast rotating disk of ro¬ 
tation speed of > 180 km/s (99% confidence level) with 
the most likely rotation velocity being 225km/s. This 
thus seems to provide an additional piece of evidence that 
damped Lya systems may be the progenitors of present 


large disk galaxies. We adopt the observed lower bound 
on the rotation speed of 180 km/s to compare with mod¬ 
els. Figure 11 shows the density of collapsed baryons 
in systems with rotation velocity (equation 18 of C98) 
greater than 180 km/s as a function of redshift, com¬ 
puted using GPM (C98). The observations are from 
Lanzetta et al. (1995b, L95 henceforth and in Figure 
11) and Storrie-Lombardi et al. 91996, S96 henceforth) 
with Icr errorbars. The errorbars are those adopted in 
S96 for both sets of observations. A small complica¬ 
tion arises from the fact that the observations are tab¬ 
ulated only for for models with qg = 0 and go = 0.5 
[Table 1 of S96]. We compute the observed values for 
other models by linear interpolating hQg between the two 
tabulated models using weighting of the inverse of the 
Hubble constant at each redshift bin. This is justified 
since the absorption distance interval is inversely propor¬ 
tional to the Hubble constant at the redshift in question: 
AX = {1 + z)'^[Ho/H{z)] (Bahcall & Peebles 1969), where 
H{z) = 77o\/f^o(l + -l- (1 — Ho — Ao)(l -I- z)^ -I- Aq 

(Peebles 1993). It is, however, still approximate due 
to non-uniform distributions of absorption redshift paths 
within each bin and finite sizes of the redshift bins. 

Without running detailed hydrodynamic simulations we 
wish only to conservatively assume that all the baryons 
(which is proportional to the mass with a fixed ratio equal 
to the global ratio) that collapse into these systems are 
in neutral form. This is conservative in the sense that it 
is an upper bound on Hg, since some amount of hydro¬ 
gen may be in other forms such as ionized or molecular 
form. Bearing in mind that the computed quantity is an 
upper bound, we see that HCDM is inconsistent with ob¬ 
servations at a high confidence level, in agreement with 
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Fig. 11.— The gas density of collapsed baryons in systems with rotation speed of > 180km/s, as a function of redshift for six models. The 
observations are from Lanzetta et al. (1995) and Storrie-Lombardi et al. 91996) with Itr errorbars. 


the conclusion drawn by Ma et al. (1997) for a compa¬ 
rable model, and consistent with previous analyses (Sub- 
ramanian & Padmanabhan 1994; Kauffmann & Chariot 
1994; Mo & Miralda-Escude 1994; Klypin et al. 1995), 
which use an HCDM model with a somewhat larger frac¬ 
tion of mass being in hot neutrinos (flh, = 0.3). Note that 
the newer observation with significantly absorption path 
length in S96 than in L95 has revised and significantly 
lowered the point at z = 3.19. But the basic conclusions 
previously reached by other authors with regard to HCDM 
model remain valid. All the other models are consistent 
with observations with perhaps the exception of OCDM25 
(open Hq = 0.25), which appears to have overproduced 
the neutral hydrogen gas in these systems. But we re¬ 
call that photoionization and other heating sources could 
significantly reduce the neutral hydrogen fraction in these 
systems. Therefore, OCDM25 cannot be ruled out yet 
without more detailed calculations. 

3.7. -High Redshift Galaxies 

Recently, observations of high redshift galaxies have 
made major break-throughs by using an innovative tech¬ 
nique to systematically and efficiently search for them 
(Pettini et al. 1997; Steidel et al. 1998, S98 henceforth). 
Interesting discoveries have already been made with a rel¬ 
atively small sample of such galaxies (~ 200). S98 find a 
large concentration of galaxies in redshift space near red¬ 
shift three in a field of size about 9 x 15 arcmin^. The 
concentration is so large for z ~ 3 that S98 are able to 
draw some interesting implications for cosmological mod¬ 
els and galaxy formation. We will use a somewhat dif¬ 
ferent method to conservatively assess the implications of 
such high redshift galaxy observations for models. In par¬ 
ticular, we will try to address how likely such a large con¬ 


centration of galaxies will occur in various models. 

We first note that Pettini et al. (1997) give a selection 
function (in redshift) of their observations simply by fit¬ 
ting the observed distribution of galaxies in the redshift 
range z = 2 — 4. This approach would be valid, if the dis¬ 
tribution of galaxies in redshift in the range of interest is 
constant. In general, the density of galaxies needs not be 
constant in redshift. Without a priori determined selec¬ 
tion function, we simply use a Gaussian selection function, 
exp (z — zYI2a'gp, where the mean z is chosen to be 3.07 
to maximize the selection function value at the peak of 
the observed galaxy concentration (z = 3.05 — 3.09) (this 
is conservative in the sense that it maximize a model’s abil¬ 
ity to match such a large concentration), and usf = 0.35 
will be used. Such a choice of selection function seems ap¬ 
propriate for the observational technique used which opti¬ 
mizes detection at z ~ 3 (Pettini 1997, private communica¬ 
tions), and should be adequate for our purpose in the sense 
that our results do not sensitively depend on the choices of 
z and cr within reasonable ranges. By the way, this chosen 
selection window function turns out to be quite close to 
what is empirically obtained from observations (Figure 2 
of Pettini et al. 1997). 

Figure 12 shows the distribution of galaxies (halos) for 
each model without (dotted curve) and with (solid curve) 
the selection function. For each model we allow ourselves 
the liberty to adjust the threshold of the virial mass of the 
halos until the number of galaxies in the redshift range 
z = 2 — 4 matches what is observed. S98 find that 180 
galaxies satisfy their color selection criteria and about 
5% of them turn out to be stars in the field SSA22 of 
size 8'.74 x 17'.64. We will simply assume that there are 
181 X (1 — 0.05) = 172 observed galaxies in the range 
z = 2 — 4, to which our models are normalized. We 
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Figure 12 



Fig. 12.— The number of galaxies (halos) per steradian per unit redshift without (dashed curves) and with (solid curves) the observational 
selection function (see text for the definition of the observational selection function). The number of galaxies (with the observational selection 
function) in each model is normalized to what is observed in the redshift range 2 = 2 — 4 by Steidel et al. (1997). 


use GPM (C98) to compute the abundance of such ha¬ 
los. Normalizing the number density of galaxies (halos) in 
each model requires to impose quite different velocity dis¬ 
persion threshold (the velocity dispersion-mass relation is 
governed by equation 18 and the virial radius-mass relation 
by equation 12 of C98) for the virialized halos as tabulated 
in Table 3. Preliminary observational evidence indicates 
that the high redshift galaxies have velocity dispersions in 
the range (T|| = 180 — 320km/s (Steidel et al. 1996). We 
thus see that the current observed range in (j| | has already 
become interesting in terms of imposing a constraint on 
the models. If future observations can significantly nar¬ 
row the range, then the combination of (T|[ and density 
may be quite powerful to discriminate between models. 
We further note the volume-limited redshift distributions 
of galaxies (solid curves in Figure 12) are quite different 
in different models. So a direct observation of this (with¬ 
out narrow observational selection window functions) in a 
sufficiently large range in redshift, perhaps by a combina¬ 
tion of several observational techniques/passbands, will be 
very important. 

The field of view of size 8'.74 x 17'.64 corresponds to dif¬ 
ferent sizes (Ta;, Ly) in comoving length units for different 
models, so does the redshift bin size (Az = 0.04; L^), all of 
which are listed in Table 3 along with some other relevant 
quantities. We compute the rms fluctuations of total mass 
in a box of comoving size Lx y- Ly x L^ for each model 
at z = 3, am,box, using the following formula in Fourier 
space: 

^ ^OO oOO pOO 

J[sin{0.5Lxkx)/i0.5Lxkx)f 

[Sm{0.5Lyky)/{0.5Lyky)]^ 

[sm{0.5L^kz)/{0.5Lzkz)f 


Pkdkxdkydkz, ( 11 ) 


where Pk is the power spectrum at z = 3. Note that we 
have assumed that there is no peculiar velocity distortion 
along the line of sight, when computing L^. This assump¬ 
tion is conservative in the sense that the actual extent in 
real space would be larger than the indicated Lx, if the 
linear compression along the line of sight for an overdense 
region is more important than the finger-of-god effect due 
to the internal velocity dispersion of the collapsed object. 
We have the following reason to believe that this is the 
case. Given the observational evidence that such peaks 
are fairly common (S98), it is highly unlikely that they 
are very rich clusters of galaxies, thus internal velocity dis¬ 
persion is unlikely to be important on the indicated scales 
of Lz- To be somewhat more quantitative, a velocity dis¬ 
persion of lOOOkm/s corresponds to 5/i“^Mpc comoving 
at z = 3 in an flo = 1 model, significantly smaller than 
Lz = 15/i“^Mpc in the same model. 

Let us now compute the bias factor of the halos at z = 3 
using the method similar to that used to compute the cor¬ 
relation of clusters (equation 28 of G98): bhaio = + 1) 

which is listed in Table 3. < f' > is computed using equa¬ 
tion 29 of G98 and oo is the rms fluctuation, both com¬ 
puted for a density held at z = 3 smoothed by a Gaussian 
window chosen to produce the correct number density of 
halos (G98). We see that peaks are moderately biased over 
mass in low density models but strongly biased in = 1 
models, in agreement with the conclusion drawn by S98. 

If a linear bias of halos over mass on the relevant large 
scales (i.e., on the box scale) exists, as supported by 
detailed galaxy formation simulations (Gen & Ostriker 
1992b; Blanton et al. 1998) and analytical arguments (Gra- 
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Table 3 

High Redshift Galaxies in Six Models 


Model 

{h ^Mpc) 

Ly [h ^Mpc) 

Lz [h ^Mpc) 

(T|| (km/s) 

^m,box 

^halo 

tCDM 

7.63 

15.4 

15.0 

165 

0.138 

4.82 

HCDM 

7.63 

15.4 

15.0 

115 

0.131 

6.38 

OCDM25 

10.8 

21.8 

22.7 

243 

0.320 

1.90 

OCDM40 

9.65 

19.5 

20.2 

224 

0.239 

2.71 

ACDM25 

11.9 

24.1 

29.3 

202 

0.207 

3.30 

ACDM40 

10.4 

20.9 

23.4 

185 

0.178 

3.92 
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Fig. 13. — Peak height of the box (yhalo box) against the bias of galaxies over the matter on the box {b^alo) is plotted for each model. The 
range in vertical axis for each model is (2cr) statistical fluctuations due to small number of galaxies in the observed galaxy peak. Also shown 
above each vertical errorbar is the velocity dispersion threshold of halos in each model, along with the observed range shown near the bottom 
of the plot. 
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mann & Einasto 1992), one obtains 


^halo,box — ^halo^m^box^ ( 12 ) 

where (Thaio,box is the rms fluctuations of halos in a box of 
comoving size Lxy< LyX Lz ai z = 3. We hrst compute the 
number of galaxies (without the selection function) in each 
model in the redshift range 2 = 2.4 — 3.8 consisting of 35 
redshift bins each of size A 2 ; = 0.04 (Steidel et al. 1998), 
fV 2 , 4 - 3 , 8 - Dividing A 2 , 4 - 3,8 by 35 gives the average num¬ 
ber of galaxies in each redshift bin of size, N. Since the 
total number of spectroscopically identified galaxies is only 
78 in the redshift range z = 2A — 3.8, we need to compen¬ 
sate for the fraction of spectroscopically yet unidentified 
galaxies, which will be assumed here to have the same red¬ 
shift distribution as the spectroscopically identified ones. 
This is accounted for by multiplying the relevant galaxy 
number by the ratio, N2a-3.s/7S. Now, we can compute 
the height of the observed galaxy peak at z = 3.07 in terms 
of the rms fluctuations in the indicated box {L^ x Ly x Lz) 
in each model: 


^box 


Npeak{N2A-3.8/78)/N 

^halo,box 


- 1 , 


(13) 


where Np^ak is the number of galaxies in the observed 
redshift peak near z = 3. There are 15 spectroscopically 
identified galaxies in the peak (Az = 3.05 — 3.09) in the 
SSA22 held, which gives a 2a upper and lower range of 
Npeak = 15 ± 8 due to Poisson huctuations. 

Figure 13 shows Vbox as a function bhaio- The range 
in the vertical axis for each model is due to {2a) Poisson 
statistical huctuations of the small number of galaxies in 
the observed galaxy peak. Also shown above each verti¬ 
cal errorbar is the velocity dispersion threshold of halos 
in each model, along with the observed range shown near 
the bottom of the plot. Clearly, the vertical ranges are 
large. Hence, a larger sample of such events will be invalu¬ 
able to dramatically shrinking the size of the errorbars. If 
such galaxy concentrations turn out to be as common as 
2a events in the real world, then all the models would 
have to be re-evaluated. The two Dq = 1 models (tCDM 
and HCDM) have one additional inconsistency with ob¬ 
servations: they do not have enough massive halos to host 
galaxies with velocity dispersion as large as indicated by 
the preliminary observations. 

It should be pointed out that the GPM method used 
here to compute the halo statistics ignores merger. But 
it may be argued that, in this case, this peak count¬ 
ing method may be more accurate than direct N-body 
simulations for precisely the reason that merger is ig¬ 
nored. Overmerging of galaxies (identified with halos) in 
N-body simulations is error-prone due to lack of realis¬ 
tic physics/sufHcient resolution. In real world, although 
halos may still merge, most galaxies are likely to retain 
their separate identities (for example, a rich cluster like 
Coma contains hundreds of galaxies but might consist of 
only a handful of dark matter halos). This is, of course, 
not in conflict with some unavoidable processes occurring 
in the real world that may transform the morphologies of 
galaxies but do not produce direct merger, such as galaxy 
harassment (Moore et al. 1996). Thus, our treatment by 
ignoring merger may in fact be closer to truth than halos 


identified in direct N-body simulations in terms of cor¬ 
rectly counting galaxies. If halo merger is more important 
in objects such as the observed galaxy concentration in 
the real world, inclusion of some merger effects would ex¬ 
acerbate the problem, as merger would reduce the number 
of galaxies in high density region (as considered here) and 
reduce the bias of galaxies (in number) over mass in the 
models. Interestingly, our results, namely that the high 
redshift galaxies are highly biased over mass, are in broad 
agreement with results from N-body simulations (Bagla 
1998; Jing & Suto 1998; Governato et al. 1998). 

We note in passing that directly measuring the bias of 
galaxies over mass at high redshift will prove to be ex¬ 
tremely difhcult. The only direct measurement of the dis¬ 
tribution of underlying mass has to come from direct map¬ 
ping by gravitational leasing. But such a measurement, 
even if possible in principle, is very difficult. This situa¬ 
tion is further worsened by the fact that matter at z = 3 
does not make an optimal contribution to the lensing effect 
compared to matter at lower redshifts. There is, however, 
one possible, though somewhat less direct, measurement of 
the bias of galaxies over mass. It makes use of the fact that 
it is possible to select a population of Lya clouds whose 
clustering strength is comparable to that of the underly¬ 
ing matter. Then, a comparison between the clustering of 
such Lya clouds and the clustering of galaxies both at the 
same high redshift would yield a measurement of bias (of 
galaxies over matter). This technique has recently been 
proposed by Gen et al. (1998). But detailed computations 
of various models are needed in order to understand the 
relationship between matter and Lya clouds in different 
models, before one can make self-consistent calibration to 
infer the clustering of mass from that of Lya clouds. 

3.8. Reionization of the Universe 

To properly tackle this issue requires detailed hydrody¬ 
namic modeling with three-dimensional radiative transfer. 
Cosmological hydrodynamic simulations (Gen & Ostriker 
1993a; Ostriker & Gen 1996; Gnedin & Ostriker 1997) 
as well as analytic analyses (Haiman & Loeb 1997) have 
demonstrated that the detailed picture of reionization of 
the universe depends on many uncertain parameters, in¬ 
cluding the properties of the ionizing sources, the dumpi¬ 
ness of the baryonic density field, the mean baryonic den¬ 
sity, etc. 

We want to address the following comparatively simpler 
question: When was the universe completely reionized? 
It is possible to obtain a relatively accurate result if one 
only wants to identify the epoch when the universe gets 
completely reionized (ignoring the detailed development 
of the reionization process). In this case one just needs to 
count the number of ionizing photons. When the number 
of ionizing photons emitted per hydrogen atom in one (hy¬ 
drogen) recombination time is equal to one, the universe 
(with regard to hydrogen) is completely ionized. We now 
formulate this criterion quantitatively, assuming that the 
universe is reionized primarily by UV ionizing radiation 
from massive stars of early galaxies. 

The Hubble time as a function of redshift is (Peebles 
1993) 

tniz) = to [Oo(l -l- z)^ -l- (1 — Hq — Ao)(l -I- z)^ -I- Ag] ^ , 

(14) 
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where to = 9.78/i ^ Gyr is the present Hubble time. The 
hydrogen recombination time is 


= 2.0 

(15) 

where 1.16 comes from the fact that the primordial gas has 
a (76%, 24%) (H, He) composition, Q;(r) is the hydrogen 
recombination coefficient and C = < pi > / < pb >'^ is 

the gas clumping factor (we have ignored here the differ¬ 
ence between total baryonic matter and ionized hydrogen 
assuming that hydrogen is ionized everywhere). For the 
second equality in equation 15 we assume T = 1.5 x lO^K, 
giving a = 6.9 x 10“^^ cm^/sec (Cen 1992). Assuming 
that a fraction, euv, of the rest mass energy of collapsed 
systems is converted into ionizing photons in a short time 
scale (i.e., star formation for the relevant massive stars 
is assumed to be almost instantaneous once a galaxy has 
collapsed), we can write down the equation, equating the 
number of ionizing photons emitted per hydrogen atom in 
a recombination time to unity, as: 


euvTTipC^ dfg dz 

hvQ dz dtn 


tree — Ij 


(16) 


where dfg/dz is the fraction of total baryonic matter that 
has collapsed and condensed out to form galaxies per unit 
redshift, rup is proton mass, c is the speed of light and hizo 
(= 13.6eV) is the hydrogen ionization potential. Using 
equations (14,15) to expand equation (16), one finds 


dfg 


^ = 1.2 X IQ-^h 'e 




flbh"^ 


l + zV^^ 

6 ) 


dz ■■ VO-0125 

[Hq -f (1 — Hq — Ao)(l + z) ^ + Ao(l + z) 


+ —(1 — Hq — Ao)(l z) ^ 


31 - 3/2 


(17) 


The main uncertainties in this treatment include: the star 
formation history of the high redshift galaxy, the initial 
mass function (IMF) of the stars, and the fraction of UV 
ionizing photons that escape host galaxies and enter the 
inter-galactic medium (IGM), the dumpiness of the gas 
and complicated interplay between radiation self-shielding 
of dense, optically thick regions and absorption of UV pho¬ 
tons by these regions. An accurate treatment of this prob¬ 
lem awaits detailed galaxy formation simulations with a 
sufficiently realistic treatment of three-dimensional radia¬ 
tive transfer. The present study serves to illustrate the 
large differences between models with regard to the reion¬ 
ization epoch, which may survive detailed, more accurate 
simulations. In the present treatment all these various 
factors have been conveniently condensed into the single 
parameter ejjv^ being the fraction of the rest mass of stars 
that turn into ionizing UV photons escaping into the IGM. 
We will use euv = (3 ± 2) x 10“^, which is lower than the 
values previously used (Gnedin & Ostriker 1997; Gen & 
□striker 1998) by approximately a factor of two. This 
factor of two is intended to account for the slight over¬ 
counting of the number of ionizing photons in equation 
16, where all ionizing photons are treated to be at Lyman 
limit, in contrast to an Orion-like spectrum of a young stel¬ 
lar association (Scalo 1986). A clumping factor C = 7.0 is 


used, which may be uncertain by a factor of a few, perhaps 
in the upward direction. 

The required galaxy formation rate per unit redshift 
dfg/dz (equation 17) is shown as the cross-shaded bands 
in Figure 14, where the upper and lower limits are due to 
the uncertainty in euv- The solid curves in Figure 14 show 
the fraction of matter in collapsed (virialized) systems per 
unit redshift with velocity dispersion greater than 35km/s 
(corresponding to virial temperature of ~ lO^Kelvin), as 
a function of redshift, that are assumed to emit the ion¬ 
izing photons. Smaller systems are unlikely to be able to 
collapse due to the lack of efficient coolants at T < lO^K 
in the presence of significant photo-ionizing radiation Held 
(Efstathiou 1992; Quinn, Katz, & Efstathiou 1996) or in 
the presence of feedback processes from young galaxies 
(Gen & Ostriker 1993a); molecular hydrogen could have 
served as an efficient coolant at the lower temperatures 
but they are likely to have been destroyed long before these 
sub-galactic size objects start to form (Ostriker & Gnedin 
1996; Haiman, Rees, & Loeb 1997). 

We see that different models predict quite different 
epochs of reionization of the universe, which are tabu¬ 
lated in Table 4 for the ease of comparison. It is noted 
that the exact epoch of reionization is quite insensitive to 
small variations of euv or the degree of dumpiness of the 
medium, simply because it is during the rapidly rising part 
of the evolution of the collapse fraction when the reioniza¬ 
tion occurs. It is clear that observations to higher redshift 
would provide a powerful discriminator between models. 
Eor example, if the universe is observed to be transparent 
at z > 5, then HGDM will be inviable; if the universe is 
observed to be transparent at z > 10, the all but the open 
models will be ruled out. 

While it is beyond the scope of current paper to con¬ 
duct detailed calculations of the properties of high redshift 
objects around the reionization epoch in each model, we 
only briefly comment that near future observations by me¬ 
ter waves telescopes such as Giant Meterwave Radio Tele¬ 
scope (GMRT; Swarup 1984) and the Kilometer Square 
Array Interferometer (SKAI; Braun 1995) and the next 
generation infrared telescopes such as the Space Infrared 
Telescope Eacility (SIRTF; Werner & Eisenhardt 1988) 
and Next Generation Space Telescope (NGST) will shed 
light on the properties of objects at the reionization epoch 
such as neutral hydrogen condensates (Scott & Rees 1990; 
Subramanian & Padmanabhan 1993; Madau, Meiksin, & 
Rees 1997) and ionizing sources (Haiman & Loeb 1997; 
Miralda-Escude & Rees 1997). In addition, the differences 
in the reionization time thus in optical depth to CMB pho¬ 
ton and in angular scales would leave different imprints on 
GMB at small angular scales. 

3.9. CMB Experiments: MAP and PLANCK 

Needless to say, future GMB experiments are likely to 
provide a very rich amount of information, that can be 
used in many aspects to constrain cosmological models 
as well as understanding astrophysical processes involved. 
However, cosmological parameters determined from GMB 
experiments alone are likely to display certain degenera¬ 
cies (e.g., Efstathiou & Bond 1998). Other astronomical 
and cosmological phenomena discussed in the preceding 
sections are, in this sense, complementary to GMB obser¬ 
vations. Here we choose to focus on one parameter, which 
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Fig. 14. — The fraction of matter formed per unit redshift in collapsed (virialized) systems with velocity dispersion greater than 35km/s 
(corresponding to virial temperature of ~ lO'^Kelvin), as a function of redshift in six models. Each cross-shaded band indicates the required 
range with which enough UV photons are produced to completely reionize the entire universe. 


Table 4 

Reionizaton Epochs in Six Models 


Model 

^rei 

tCDM 

5.8 

- 8.2 

HCDM 

3.7 

- 5.6 

OCDM25 

28.0 

- 36.0 

OCDM40 

13.0 

- 20.0 

ACDM25 

10.0 

- 13.0 

ACDM40 

7.7- 

- 10.0 
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Fig. 15.— summarizes the model results on n when all the models are normalized to both COBE and local clusters, assuming that the 
primordial spectrum has a pure power law shape for the scales of cosmological interest (~ 1 — 6000h“^Mpc). The ranges in n are 2cr. Also 
shown are the 2cr errorbars that MAP and PLANCK are thought to be able to achieve (Zaldarriaga, Spergel, &; Seljak 1997). Note that the 
horizontal locations of MAP and PLANCK indicators are arbitrary. 


is cosmologically important but is not tangled with other 
complicated processes such as recombination and reion¬ 
ization processes, namely, the primordial power index n. 
Anot her reason for favoring the choice of n is that n is 


distin :tly different between models and perhaps offers one 
of cle anest tests of models. The underlying reason that 


differi nt models have distinctly different n is a product 
of combined effort: combining CMB (COBE) with local 
cluster observations. It is clearly seen from Table 1 that 
normalizing a model both to very large scale (COBE) and 
cluste r scale results in a tight constraint on the primor- 


dial p |ower spectrum index n. However, currently COBE 
does not provide strong constraint on n due to its large 
errorbars [An = 0.6 (2 (t)]. In the foreseeable future there 
are two satellite missions, MAP and PLANCK, which are 
proposed to produce much higher resolution (arminute to 
sub-degree versus ~ 10° of COBE) and much more sensi¬ 
tive CMB temperature fluctuation measurements. 

We summarize the model results on n in Figure 15. The 
ranges in n are 2a. Also shown are the 2a errorbars that 
MAP and PLANCK are thought to be able to achieve (Zal¬ 
darriaga, Spergel, & Seljak 1997). Note that the horizontal 
locations of MAP and PLANCK indicators are arbitrary. 
It is seen that MAP will be able to provide a powerful dis¬ 
criminator between the models, whereas PLANCK should 
be able to constrain n thus test models with exquisite pre¬ 
cision. 

4. DISCUSSION 

When models are normalized to both COBE on very 
large scales and clusters of galaxies on intermediate-to- 
large scales, They naturally have grossly similar properties 
of clusters of galaxies at zero redshift. Other local tests 


surrounding clusters regions or weighted heavily by clus¬ 
ters are unlikely to be very sensitive for the purpose of dif¬ 
ferentiating between models. Among suc h tests are galaxy- 


galax y pairwise velocity dispersion (e.g., [Davis fc Peeblei 


1983; Strauss, Ostriker, & Cen 1998), galaxy-ga l axy two 


point correlation function (Davis fc Geller 1976; Peeblei 


1980 ) and large-scale velocity flows in the local universe 
(e.g., Strauss et al. 1995). Although it is conceivable that 
fine details within the clusters at zero redshift may dif¬ 
fer among different models, such as cluster morph ology, 
substructur e etc. (Richstone, Loeb, & Turner 1992; Mohi 


et al. 1995; see Cen 1997b for a somewhat more realistic 


view on this considering some unavoidable projection ef¬ 
fects) , it is likely that more significant differences between 
the models will come from regions far removed from local 
clusters of galaxies. There are two ways to be distant from 
clusters. First, one can move to regions such filaments 
or voids in the local universe which are spatially distant 
from the clusters. There is no doubt that there are signif¬ 
icant differences in the properties of these regions among 
different models. But detailed modeling of such regions, 
in particular, voids, is extremely difficult due to the lack 
of appropriate small scale resolution in such regions even 
with the best simulations available. We have followed on 
the second route in this paper by traveling back in time 
away from the local universe and studied the properties 
of model universes at moderate-to-high redshift through a 
suite of cosmological phenomena. 

In addition to a consideration of galaxy power spectrum 
from SDSS galaxy redshift survey (§3.2), which is shown to 
be capable of differentiating between models, we have ex¬ 
amined eight phenomena ranging in redshift from zero to 
that of last scattering. They include the correlation func- 
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tion of very rich clusters of galaxies (§3.1), the evolution of 
the abundance of rich clusters of galaxies (§3.3), the grav¬ 
itational lensing by moderate redshift clusters of galaxies 
(§3.4), the properties of Lya forest at high redshift (§3.5), 
the damped Lya systems at high redshift (§3.6), the high 
redshift galaxies (§3.7), the reionization epoch of the uni¬ 
verse (§3.8) and the future high resolution, high sensitivity 
CMB experiments (§3.9). Some of these observations have 
already provided discriminating tests between some mod¬ 
els. For example, HCDM model is strongly ruled out by 
the observed large density of neutral hydrogens in damped 
Lya systems and possibly inconsistent with observations 
of Lya clouds. Both Oq = 1 tCDM and HCDM mod¬ 
els are in disagreement with the observed abundance of 
rich clusters of galaxies at z ~ 0.4 — 0.8 and the observed 
strong clustering of high redshift galaxies at z ~ 3, while 
other, low density models are in fair agreement with those 
observations. 

We find that all these moderate-to-high redshift phe¬ 
nomena possess potentially very powerful capabilities of 
differentiating between all the remaining “viable” mod¬ 
els. We now discuss the requisite condition for each phe¬ 
nomenon to be capable of providing some useful con¬ 
straint. 

It is found that different models differ most for corre¬ 
lations of very rich clusters of galaxies, those with mean 
separation > 90/i“^Mpc. Let us take the data point shown 
as a solid dot in Figure 3 as an example to quantify this. 
If the errorbar were 2/i“^Mpc (Icr) (with the mean value 
unchanged), all the models may be ruled out at 4(j confi¬ 
dence level. Is this achievable? The answer is a probable 
yes with SDSS. Let us assume that the correlation function 
is a powerlaw with a slope of —1.8. Then, the mean corre¬ 
lation function has a value ^(ro) = 3^(ro)/(3 — 1.8) = 2.5 
at ro (where the correlation function has a value unity). 
In order to achieve a statistical errorbar of 2/i“^Mpc for 
ro = 42/i“^Mpc, Npair = (42/2)^ = 442 pairs of clus¬ 
ters are necessary. Since f = Npair/(Vsurveyn) — 1 = 2.5 
at ro, we get Vsurvey = Npair/3.5n = 1.0 x 10®/i“^Mpc^ 
[n = (l/94)3/i3Mpc-3 is assumed]. The SDSS survey vol¬ 
ume is about ^500^ = 5.2 x 10®ft,“^Mpc^. So it seems 
that clusters from SDSS will provide an unprecedented 
accuracy for the determination of cluster correlations of 
very rich clusters. 

The evolution rates of the abundances of rich clusters 
of galaxies in different models are quite different and dif¬ 
ferences of the abundances of clusters increasingly devi¬ 
ate from one another towards high redshift (Figure 7). 
As shown earlier, even the currently available very small 
data set is already able to differentiate the flo models 
(tCDM,HCDM) from the rest. However, larger cluster 
samples will be required to constrain Oq to a high ac¬ 
curacy. As shown §3.3, a high redshift X-ray cluster sur¬ 
vey covering a quarter of the sky with enough sensitiv¬ 
ity to measure the temperature of hot luminous clusters 
(fcT > 7keV) at z > 1 will be able to determine flo with 
an uncertainty of AOq = 0.02 (Icr), if Dq -I- Aq = 1, and 
Aflo = 0.008, if Aq = 0. These accuracies are comparable 
to those which will be determined from the next genera¬ 
tion microwave background fluctuation experiments (e.g., 
Zaldarriaga, Spergel, & Seljak 1997), thus provides a very 
powerful cosmological test and also a very important con¬ 


sistency cross-check. 

The gravitational lensing optical depths for large split¬ 
tings or giant arcs due to moderate-to-high redshift clus¬ 
ters of galaxies increasingly differ towards high redshift be¬ 
tween models. While this is closely related to the evolution 
of abundance of rich clusters, the differences between dif¬ 
ferent models tend to be larger and more extreme because 
of the highly nonlinear dependence of gravitational lensing 
cross section on cluster mass. For background sources at 
z < 1, the primary differences occur between models with 
different flo! models with a same Hq but with and without 
a cosmological constant give comparable results (Figure 
9a). For background sources at z ~ 3, all the models be¬ 
come non-degenerate (Figure 9b). A systematic search of 
about 10® background sources at z > 3 will be desirable. 
In fact, SDSS should provide such a sample (~ 10® quasars 
mostly at high redshift with photometric redshifts and 10® 
quasars with spectroscopic redshifts, Margon 1998). One 
can use such a sample in two ways. Firstly, the total num¬ 
ber of lensed sources will be different between the mod¬ 
els. With 10® quasars (assuming to be at z^ = 3 for this 
simple illustration), one should see (1900 ± 44, 710 ± 27, 
570±24, 330 ±18, 100 ±10, 39 ±6) multiple-image quasars 
with separation > 30” in (OCDM25, OCDM40, ACDM25, 
ACDM40, tCDM, HCDM), respectively. This would al¬ 
low an unambiguous test of these six models. It would 
put a constraint on Oq with AOq ~ 0.01 if the universe 
is open or AOq ~ 0.03 if the universe is flat. Secondly, 
and more sensitively, the distributions of lenses themselves 
in different models will be very different. For example, 
one expect to see (200.0, 4.0, 0.5, 0.005, 0.0, 0.0) lenses 
at z > 2 in (OCDM25, OCDM40, ACDM25, ACDM40, 
tCDM, HCDM), respectively, that produce arcs or split¬ 
tings > 30”. 

The properties of Lya forest at high redshift are likely 
to be dramatically different among the models given the 
order of magnitude different powers on the relevant scales. 
With regard to this, observations are perhaps slightly 
ahead of theory; there is a large and still rapidly accu¬ 
mulating observational database for Lya forest. But de¬ 
tailed hydrodynamic simulations of Lya forest, while hav¬ 
ing made ground-breaking advances in our understanding 
of Lya forest phenomenon in recent years, are yet to sort 
out the fine dependence of results on such issues as vari¬ 
ations among different hydrocodes, resolution and boxsize 
effect, radiative transfer effect and star formation feedback 
effects, before they can be utilized to test cosmological 
models in accurate (thus meaningful) ways. Preliminary 
comparisons between models and observations using the 
one-point flux distribution probability function (Rauch et 
al. 1997) favor the two A based models (ACDM25 and 
ACDM40), but disfavor HCDM and OCDM25 in the oppo¬ 
site sense that HCDM has too little power and OCDM25 
has too much power on the relevant scales. tCDM and 
OCDM40 appear to be marginally consistent. 

The neutral hydrogen mass in damped Lya systems at 
high redshift is undoubtedly a very sensitive test. Their 
role is comparable to that of rich clusters today in the 
sense that damped Lya systems are rare systems at the 
high redshift. Four improvements will dramatically en¬ 
hance its discriminatory power. First, more detailed kine¬ 
matic studies of these systems, at the highest redshift pos¬ 
sible, are indispensable to providing more tight limits on 
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the mass of these systems. Second, a significant reduction 
of observational errorbars with more samples can signif¬ 
icantly increase the sensitivity of the test. SDSS should 
provide a much larger quasar sample. Third, it is beneficial 
to go to higher redshift to further increase its sensitivity of 
testing cosmological models. Finally, detailed modeling of 
these systems is necessary in order to determine the pre¬ 
cise amount of gas in neutral form. For this we will have 
to implement detailed three-dimensional radiative trans¬ 
fer in the existing hydrocodes, run large simulation boxes 
(> 20/i“^Mpc) to ensure a fair sample and have a bet¬ 
ter understanding of the star formation processes in these 
systems. 

The abundance and clustering of high redshift galaxies 
at z ~ 3 has been shown to differ among models. These are 
the large proto-clusters or proto-superclusters that even¬ 
tually produce the clusters and superclusters of galaxies 
that we see today. The issue here is to enlarge the ob¬ 
servational sample size. If the observed galaxy concentra¬ 
tion seen by S98 turns out to be ubiquitous, all the mod¬ 
els considered would be in a very uncomfortable position. 
With about ten fields like SSA22, preferably randomly dis¬ 
tributed spatially, we may be able to reject some models 
with high confidence levels. A perhaps more discrimina¬ 
tory test with the high redshift galaxies is from measuring 
the velocity dispersion (or mass) of these galaxies. If the 
high redshift galaxies are as massive as the present day 
L* galaxies, as present preliminary observational evidence 
indicates (Steidel et al. 1996), all Hq = 1 models are invi- 
able. 

The reionization epoch of the universal baryons is shown 
to be substantially different among the models. It is very 
worthwhile to search for quasars at as high redshift as pos¬ 
sible. For example, if quasars can be seen beyond z > 5, 
HCDM will be inviable; if quasars can be seen beyond 
z > 10, only low density open models could pass. Next 
generation infrared telescopes such as SIRTF and NGST 
and next generation radio telescopes such as GMRT and 
SKAI will prove to be extremely useful to not only test 
the models but also understand the detailed history of the 
reionization. Future CMB experiments will also shed light 
on this. 

Finally, the future high resolution, high sensitivity CMB 
experiments such as MAP and PLANCK will provide 
many sensitive tests. Among them one critical test will be 
an unambiguous determination of the primordial power in¬ 
dex. Accurate determinations by MAP and PLANCK will 
provide strong constraints on 

5. CONCLUSIONS 

We examine a set of observations and compare them (if 
currently available) with a set of six COBE and cluster 
normalized CDM models, or forecast their capabilities in 
constraining models or outline the requirements for them 
to be capable of constraining models in a dramatic fash¬ 
ion. Detailed discussion has been presented in the previ¬ 
ous section. In addition to SDSS galaxy redshift survey 
(§3.2), we have examined eight different phenomena in¬ 
cluding the correlation function of very rich clusters of 
galaxies (§3.1), the evolution of abundance of rich clusters 
of galaxies (§3.3), the gravitational leasing by moderate- 
to-high redshift clusters of galaxies (§3.4), the properties 
of Lya forest (§3.4), the damped Lya systems at high red¬ 


shift (§3.6), the high redshift galaxies (§3.7), the reioniza¬ 
tion epoch of the universal baryons (§3.8) and the future 
high resolution, high sensitivity CMB experiments (§3.9). 
Some of these observations have already provided discrim¬ 
inating tests on some models. 

The hot-l-cold dark matter model (with = 0.2) is 
already ruled out at a > 3cr level by current observa¬ 
tions of damped Lya systems at high redshift and the 
observed high abundance of rich clusters of galaxies at 
z ~ 0.4 — 0.8. There is some preliminary indication that 
there is not enough power in HCDM to account for the 
observed high redshift galaxies and the Lya forest lines. 
HCDM is marginally consistent with the existence of z ~ 5 
quasars in the sense that the HCDM universe would be 
marginally completely reionized by z ~ 5; seeing higher 
redshift quasars will be an embarrassment to HCDM. 
Overall, HCDM is the least viable model among exam¬ 
ined and is ruled out in several independent ways at high 
confidence levels. 

The tilted CDM model (tCDM) is in disagreement with 
the observed high abundance of rich clusters of galaxies at 
z ~ 0.4 — 0.8 at a high confidence level (> 3cr). It is in¬ 
consistent with the observed strong correlation of clusters 
of galaxies — another manifestation of lack of large-scale 
power of the model (even with a substantial tilt, as is re¬ 
quired to fit both COBE and local cluster abundance and 
included in the analysis). tCDM is marginally inconsistent 
(~ I — 2a) with the observations of damped Lya systems 
in the sense that it underpredicts the amount of neutral 
hydrogen in damped Lya systems at z ~ 3 — 4, but is 
marginally consistent with the observed flux distribution 
of Lya forest. In summary, tCDM does not seem to be 
favored. 

Low density, flat CDM models with a cosmological con¬ 
stant are consistent with all the observations examined 
here with the preferred range of the density parameter be¬ 
ing Do = 0.30 ± 0.10. This preferred range is primarily 
constrained at the high end by the abundance of clusters 
of galaxies at z ~ 0.5 (in order not to underproduce the 
cluster abundance) and at the lower end by the flux distri¬ 
bution of Lya forest at z ~ 2 (in order not to overproduce 
the fluctuations in the flux distribution) and also by ob¬ 
servations of the damped Lya systems z ~ 3 (in order not 
to overproduce the abundance of damped Lya systems). 

Low density, open CDM models are consistent with 
all the observations examined that are currently available 
with the preferred density range being Dq = 0.45 ± 0.10, 
which is constrained by similar observations aforemen¬ 
tioned in the preceding paragraph for the ACDM mod¬ 
els. An additional constraint is provided by the observed 
galaxy fluctuations in the local universe, as{gal) ~ 1.0. 
Since it is improbable for galaxies to be anti-biased trac¬ 
ers of the mass, a^{gal) ~ 1.0 (Strauss & Willick 1995) 
implies as (mass) < 1.0, which translates into an require¬ 
ment Do > 0.25 for open CDM models. 

Some improvements on the comparisons between models 
and observations will be likely in a relatively short term, ei¬ 
ther through more realistic modeling and/or by improved 
observational data, including the correlation function of 
very rich clusters of galaxies at z=0 (from SDSS), the Lya 
forest at z ~ 2 — 4 (primarily by more careful model¬ 
ing and analysis), damped Lya systems (by Keck) and 
high redshift galaxies at z ~ 3 (using Keck and others). 



24 


These observations look at different physical scales of a 
model at different redshift: the correlation function of very 
rich clusters of galaxies probes the power on large scales, 
A > 100/i“^Mpc, at z = 0; the Lya forest lines measure 
the power on small scales, A ^ 0.1/i“^Mpc, at z ^ 2 — 4; 
the damped Lya systems check the power on intermediate 
scales, A ~ 1.0h“^Mpc, at z ^ 3; and the abundance and 
clusters of the high redshift galaxies constrain the power 
on intermediate-to-large scales, A ~ 1.0 — 100/i“^Mpc, at 
z ~ 3. These improvements will further tighten the pa¬ 
rameter space for the low density models. 

The evolution rates of the abundances of rich clusters 
of galaxies in different models are quite different and dif¬ 
ferences of the abundances of clusters increasingly devi¬ 
ate from one another towards high redshift. A large and 
deep X-ray cluster survey would be invaluable to deter¬ 
mine the abundance of high redshift, hot and luminous 
X-ray clusters, which would permit an accurate determi¬ 
nation of flo- A large survey covering about a quarter of 
the sky and capable of detecting hot and luminous clus¬ 
ters [kT > 7keV) at z > 1.0 would place a constraint on 
Ho with AHo ~ 0.01. These accuracies are comparable 
to those which will be determined from the next genera¬ 
tion microwave background fluctuation experiments (e.g., 
Zaldarriaga, Spergel, & Seljak 1997), thus provide a very 
powerful test and also a very important consistency cross¬ 
check. While AXAF will provide a deep sample of X-ray 
clusters of galaxies, a larger survey is required to make a 
competitively accurate constraint on Hq. 

Gravitational lensing observations of giant arcs or very 
large splitting multiple quasars by clusters of galaxies will 
provide very sensitive tests. For the same reason that 
cluster abundances at high redshift are different in dif¬ 
ferent models, gravitational lensing provides an indepen¬ 
dent measurement of the same objects - rich clusters of 


galaxies at moderate to high redshift - as well as lensed 
high redshift sources. The quasar data set from SDSS may 
permit an unprecedented, accurate determination of both 
the lensing optical depth and the redshift distribution of 
the lenses themselves (perhaps in conjunction with follow¬ 
up observations to determine the redshifts of high redshift 
lenses). A very tight constraint can be placed on Hq with 
AHo ~ 0.01. 

Observation of the reionization of the universe would 
provide the potential to differentiate between open and A 
CDM models. Perhaps the simplest discriminator will be 
to sight a very high redshift quasar z > 10. If such a 
quasar is seen, open models will be preferred. 

A combination of several (or all) of the aforementioned 
observations will offer several ways to break the current 
degeneracy of models that are both COBE and cluster 
normalized, primarily between open CDM and A CDM 
models. On the other hand, any conflict among the obser¬ 
vations with respect to all the models may have profound 
implications with respect to the Gaussian picture of struc¬ 
ture formation and/or the general gravitational instability 
theory. 
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